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ABSTRACT 

The central engine causing the production of jets in radio sources may work intermit- 
tently, accelerating shells of plasma with different mass, energy and velocity. Faster 
but later shells can then catch up slower earlier ones. In the resulting collisions shocks 
develop, converting some of the ordered bulk kinetic energy into magnetic field and 
random energy of the electrons which then radiate. We propose that this internal shock 
scenario, which is the scenario generally thought to explain the observed gamma-ray 
burst radiation, can work also for radio sources in general, and for blazar in particular. 
We investigate in detail this idea, simulating the birth, propagation and collision of 
shells, calculating the spectrum produced in each collision, and summing the locally 
produced spectra from those regions of the jet which are simultaneously active in 
the observer's frame. We can thus construct snapshots of the overall spectral energy 
distribution as well as time dependent spectra and light curves. This allows us to 
characterize the predicted variability at any frequency, study correlations among the 
emission at different frequencies, specify the contribution of each region of the jet to 
the total emission, find correlations between flares at high energies and the birth of 
superluminal radio knots and/or radio flares. The model has been applied to quali- 
tatively reproduce the observed properties of 3C 279. Global agreement in terms of 
both spectra and temporal evolution is found. In a forthcoming work, we explore the 
constraints which this scenario sets on the initial conditions of the plasma injected in 
the jet and the shock dissipation for different classes of blazars. 

Key words: gamma rays: bursts — X-rays: general — galaxies: active 



1 INTRODUCTION 

Recent developments in the blazar field, mainly stimulated 
by the discovery that these sources are strong 7-ray emit- 
ters, are offering new and important clues for understand- 
ing relativistic jets. All blazar spectral energy distributions 
( SEP) are charact erized by two broad emission peaks e.g. 
(Fossati et al. 1998), strongly variable on different timescales 
(Wagner & Witzel 1995; Maraschi et al. 1994; Hartman et 
al. 1996). The location of these peaks (their frequency, rel- 
ative and absolute flux) and their variability behavior are 
a strong diagnostic tools to discriminate among theoreti- 
cal models and find the intrinsic physical parameters of the 
source. The broad band spectrum from the far infrared up- 
wards can be explained by a "single-zone" model, where the 
two peaks are respectively due to synchrotron and inverse 
Compton scattering from the same relativistic electron pop- 
ulation (Maraschi, Ghisellini & Celotti 1992; Sikora, Begel- 
man & Rees 1994) . Since the radiative cooling time is much 



shorter than the dynamical one, the emitting particles must 
be accelerated in situ all along the jet, by means of a process 
not yet univocally determined. 

The high energy fiux varies generally on a timescale of 
the order of a day, and the b ulk Lorentz factor (B LF) F of 
the jet is estimated around 10 ( Ghisellini et al. 1993| ), leading 
to an upper limit for the source dimension of the order of 
R ~ 10^^ cm. Furthermore, the required transparency for 7- 
rays to the process of photon-photon pair production, sets 
a lower limit in the most powerful sources, of order « 10^^ 
cm, to the typical distance where the bulk of the radiative 
dissipation occurs. 

Blazars are also characterized by strong radio emission, 
with a flat spectrum [F{u) oc v'" , with a between ~ —0.5 
and ~ 0.5] and variability time scales of the order of weeks, 
months (see e.g. Ulrich, Maraschi & Urry, 1997). While the 
emission at frequencies from the far infrared upwards can 
be consistently due to a single population of electrons in 
the inner parts of the jet, the emission in the radio band 
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has to be generated - in order not to be self-absorbed - 
on the outer, parsec scales where the plasma properties are 
expected to be different. The bulk motion is still relativistic 
at such distances, with BL Fs similar to those i nferred on the 
sub-parsec scales, F ~ 10 (Zensus et al. 1997). 

The existing models, both homogeneous (e.g. Sikora, 
Begelman & Rees 1994; Inoue & Takahara 1996) and in- 
homogeneous ones (e.g. Celotti, Maraschi & Treves 1991; 
Ghisellini & Maraschi 1989), estimate the jet plasma prop- 
erties, the magnetic field B and the random Lorentz fac- 
tors (RLE) of the electrons emitting at the peak of the 
spectra, 7p, by fitting the observed SED. In the inner jet 
B ~ 1 G and 7p varies from ~ 10^ to lO'' depending on the 



power of the source (Ghisellini et al. 199i): for the powerful 
OVV quasars - with bolometric luminosity of the order of 
L « 10*'' — 10"** erg s~^ - the synchrotron peak is at infrared- 
optical frequencies and 7p is estimated around 10^~^, while 
for the less powerful sources, like the weak BL Lacs - with 
L ~ 10'*'' — IC'^erg s~' - the synchrotron emission peaks 
in the UV-X-ray range and the RLFs of the electrons must 
be higher, 7p ~ 10''"''. This trend of decreasing peak en- 
ergy with increasing source power appears to be systematic 
and can be accounted for by a corresponding increase in 
the intensity of a photon field external to the jet, which 
enhances the inverse Compton cooling losses (Fossati et al. 
1998; Ghisellini et al. 1998). In the parsec scale jet the ob- 
served (radio) emission is believed to be synchrotron and 
the corresponding magnetic field and electrons RLE are es- 
timated to be i? « 10~^ G and 7p « 100, leading to particle 
cooling timescales much longer than the dynamical one. 

Although these studies have allowed us to consistently 
derive the physical parameters in the emitting region, many 
aspects in the understanding of relativistic jets remain open, 
most notably the jet energetics and particle acceleration. In 
order to explore these issues and their relationship we con- 
sider a scenario where the plasma characteristics are not 
treated as free parameters, but are the result of the jet dy- 
namics, thus relating the observed emission properties to the 
transport of energy along the jet. We achieve that by quan- 
titatively considering a general scenario for the interpreta- 
tion of blazars analogous to the standard scenario proposed 
to explain gamma-ray bursts, which was actually originally 
proposed for blazars more than two decades ago (Rees 1978) . 

The most important assumption in this scenario is that 
the central power engine produces energy which is channeled 
into jets in an intermittent way, though such a time depen- 
dent process cannot be easily inferred from first principles. 

If the physical conditions determining the energy depo- 
sition are not steady, the distribution of BLFs and masses 
is non-uniform within the ejecta: faster portions of the fiow 
would then catch up with slower ones leading to the de- 
velopment of relativistic shocks. The shocks would plausi- 
bly heat the expanding ejecta, generate/amplify a tangled 
magnetic field and accelerate leptons to relativistic ener- 
gies, in turn causing synchrotron and inverse Compton emis- 
sion. In this scenario - named internal shock model (e.g. 
Rees & Meszaros 1994; Lazzati, Ghisellini & Celotti 1999, 
Panaitescu, Spada & Meszaros 1999) - the following obser- 
vational characteristics can be naturally explained: 

• Efficiency — The radiative output of radio sources in 
general, and blazars in particular, has to be a small frac- 



tion of the energy transported by the jet (less than 10%): 
extended radio structures require in fact a power input ex- 
ceeding what is emitted (e.g. Rawlings & Saunders 1991; 
Celotti & Fabian 1993). Low efficiency is indeed a charac- 
teristic of the internal shock model, as the range of the BLFs 
of the colliding shells can be only small (say the ratio of the 
maximum and minimum ones FM/Em ~ 3) in order to have 
a resulting jet bulk motion on larger scale with E ~ 10. 

• Minimum distance for dissipation — As mentioned 
in the Introduction, there is a minimum dimension for the 
7-ray source and the bulk of the observed power must be 
produced at some distance from the jet apex, i.e. from the 
accretion disk. In other words, the emitting zone must be far 
away from X-rays sources and not extremely compact itself 
otherwise, quite inevitably, the high energy 7-rays are ab- 
sorbed by the dense radiation fields. The electron-positron 
pairs generated by such photon-photon interaction would 
then re-process the 7-ray power into softer (especially X- 
ray) radiation, which is not observed (Ghisellini & Madau 
1996). In the internal shock model the jet becomes radia- 
tive on scales much larger than the central source dimension 
{Ro ^ 10'''cm for a black hole of M = 10^ Mq) since aU the 
velocities are close to c and therefore the shells propagate 
for a significant distance before collisions can start. The ex- 
pected minimum distance for the 7-rays emission is of the 



order of R ^ RoT'' 



10' 



10 cm, where Ro/c measures 



the time interval between two consecutive shell ejections. 

• Variability — Clearly the observed large amplitude 
variability must be explained by a non steady-state model: 
we believe that the shell-shell collision scenario studied here 
is the simplest among them. 

• Continuous energy deposition along the jet — Dis- 
sipation is observed to occur along the whole jet, from the 
inner to the parsec scales with decreasing luminosity and 
typical photon frequency. The dynamical evolution of a wind 
in the internal shock model can indeed reproduce such be- 
havior: shell-shell collisions take place along the whole jet, 
leading to a continuous conversion of bulk energy into inter- 
nal one over a large range of distances. Due to the averaging 
of the shell BLFs during the expansion, the collisions be- 
come progressively less efficient as the distance increases, 
thus progressively reducing the released internal energy, the 
electron RLE, and the magnetic field. 

As this scenario allows to predict a general and fairly 
complete description of the entire radiative jet, we believe 
that it can be a powerful tool to understand blazar spectra, 
their variability and the possible correlations between ob- 
served quantities on different scales. While we plan to fully 
explore the large potentialities of this model in a forthcom- 
ing work, in this paper we present the basic assumptions and 
the general set up of the model. For illustration purposes 
(i.e. not aiming to determine precise values of the physical 
quantities) we also consider the case of a specific well stud- 
ied blazar, specifically 3C 279. Section 2 therefore describes 
the dynamics of the wind and the shocks forming from shell 
collisions, while in Section 3 the radiative properties of the 
shocked plasma are presented. The numerical simulations for 
the time dependent evolution of the system are described in 
Section 4, while the results and the specific case of 3C 279 
are the subject of Section 5. In the final Section 6 we discuss 
our findings. 
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2 HIERARCHICAL INTERNAL SHOCKS 

The entire jet emission is simulated by adding pulses radi- 
ated in a series of internal shocks that occur in a unsta- 
ble relativistic wind (Daigne & Mochkovitch 1998; Spada, 
Panaitescu & Meszaros 2000). After (parametrically) set- 
ting the dynamics of the wind ejection, we calculate the 
radii where the shells collide. For each collision we study 
the hydrodynamics, in order to determine the shock veloc- 
ity, the compression ratio and the internal energy E^h of 
the shocked fluid. Assuming a given partition of E^h among 
protons, electrons and magnetic field we calculate the rel- 
evant physical parameters in the shocked fluid, and finally 
the features of the emitted radiation. 



2.1 Ejection features 

The wind is discretized as a sequence oi N = t^/tv shells, 
where is the duration time of the wind ejection from the 
central source and <^ is the average interval between 
consecutive ejections. Each shell is characterized by a mass 
Mj, a bulk Lorenz factor Vj and an ejection time tf. at 
t = tj the j-th shell is at a distance Rj = i?o, where Ro is 
the central source dimension (i?o = lO^^cm). 

The shell BLFs Tj are chosen randomly from a uniform 
distribution between Fm and Fm, and the shell masses are 
also randomly picked up from a uniform distribution with 
an average value Mt/N, where the total mass in the wind 
Mt is determined by the requirement: 



j=i 



(1) 



where L„ is the wind luminosity. The time interval between 
the ejection of a couple of consecutive shells is assumed to 
mimic a scenario in which the inner engine behaves as an 
accumulation mechanism: the quiescent time between the 
ejection of the j*'' and the (j + 1)*'' shells is proportional to 
the total energy in the {j + 1)"" sheU. 

The ejected shell width Aj is of the same order of the 
central source dimension Ro and the shell internal energy Uj 
is negligible (all the energy is stored as bulk kinetic energy). 



2.2 Wind dynamics 

The dynamics of the wind expansion is characterized by a 
series of two-shell collisions in which the faster shells catch 
up with the slower ones in the outer parts of the ejecta. 
The calculation of the collision radii follows the procedure 
of Kobayashi, Piran & Sari (1998). Given the wind config- 
uration after the i-th collision, i.e. shell radius R^{j), shell 
velocity /3'(j), shell width A'(j) and expansion velocity of 
the shell I3l{j), we calculate the minimum value among all 
the collision times between pairs of shells {j,j — 1) with 



5ti+i — min(tjj_i) 

R{j) - R{j - 1) - 0.5(A(j) + A{j - 1)) 



1) 



/3(j)+0.5(/?e(j)+/3c(j-l))' 



(2) 
(3) 



The radius of the (i+l)-th collision represents the position of 
the two shells colliding after 5ti+i. This procedure is iterated 



for a wind described by an up-dated configuration, with a 
merged shell instead of the two colliding ones. 

Regarding the two-shell interaction as an inelastic colli- 
sion, the momentum and energy conservation laws yield the 
physical properties of the merged shell. If (F^, Mi, rji) and 
(Fo, Mo, ?7o) are the BLFs, rest masses and internal energies 
of the fast (inner) and slow (outer) shell with Fi > Fo, the 
merged shell has a rest mass Mm = Mi + Mo, a BLF given 
by: 



^iVi + poFo 

Mi/ri + Po/Fo 



1/2 



(4) 



where /li — Mi -\- rii/c? and /lo = Mo 4- rjo/c? ■ The total 
internal energy of the merged shell will then be: 

Sin = + ??o + lJi^c{V, - FnO -|- MoC^(Fo " F^). (5) 

The dynamical efficiency, namely the newly generated 
internal energy over the total energy, is given by: 



Cd = 



Mi (Ft — Fm) + Mo (ro — Fm) 
Fi Mi + Fo Mo 



(6) 



This efficiency decreases with the ratio Fi/Fo and is maxi- 
mized for Mo = M»- 

The width of the shells spreads between two consecutive 
collisions with a velocity: 



1 



(7) 



F2 1 - [PPLY ' 

as measured in the observer's frame, where pi^c — v'^ is the 
sound velocity of the plasma in the shell co-moving frame: 



1 



3 A/si 



(8) 



During the free expansion between two subsequent collisions, 
each shell loses internal energy for adiabatic losses, while the 
BLF is held constant. The magnetic field, on the contrary, is 
generated at each collision without keeping memory of the 
field generated in precedent collisions (if any). The plasma 
thermal energy Eth increases in each collision - by the frac- 
tion of Ein which is not used to either accelerate electrons 
to a power-law energy distribution and/or generate/amplify 
magnetic field - leading to an increase in the plasma sound 
velocity vi, during the wind expansion. 

Since during each collision the propagation of the shocks 
through the two shells compresses the width of the merged 
shell Am, the calculation of Am itself requires a hydrody- 
namic treatment of the interaction, that is presented in the 
next section. 



2.3 Hydrodynamic of the two— shell interaction 

Let us consider the dynamics of the interaction between two 
shells. Here it will be assumed that the presence any mag- 
netic field is negligible from the dynamics point of view. 
During the collision two shocks are formed, separated by a 
contact discontinuity: a forward shock (FS) that propagates 
in the external slow shell and a reverse shock (RS) that 
propagates backward in the faster shell (see Figure 0). 

The calculation of the shocked plasma BLF, F^h, in the 
l ab frame can be done apply ing the following conditions 
( Panaitescu fc Meszaros 199£| ): 
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Figure 1. Diagram of a two-siiell collision. A inner fast shell with 
BLF Ti and density pi catches up with a slower shell characterized 
by BLF To and density po- The reverse shock (RS), propagating 
in the inner shell, and the forward shock (FS), propagating in 
the outer shell, are shown. Both of the shocked fluids move with 
BLF Fsh in the lab frame. In the un-shocked fluid frame the RS 
BLF is Fsii i and the FS BLF is Fg^.o- The co-moving pressures 
of the shocked fluids are equal around the contact discontinuity, 
(Psh)i = (Pah)o. 



• The Blandford-McKee (1976) jump equations for the 
FS and RS. For the FS the energy and particle number con- 
servation yield: 



— ) = mpC^(rsh,o — 1) 



TFsh.o + 1 
7-1 ' 



(9) 



(10) 



where ggh and are the internal energy and the numerical 
baryon density of the shocked plasma, respectively, and n is 
the pre-shocked numerical density. All these quantities are 
measured in the fluid co-moving frame. 7 is the adiabatic in- 
dex and Fsh.o is the BLF of the shocked plasma in the frame 
of the outer shell. A similar set of equations holds for the 
RS. The shocked plasma pressure p = (7 — 1)m is calculated 
through Eqs. (||) and (|l^): (psh)o{i) = (rsh,o(i)-l)(7rsh,o(i) + 
l)Po(i)C^, with Po(i) being the co-moving matter density of 
the shells before the collision. 

• The condition of equal pressure for the shocked fluids 
around the contact discontinuity (CD) (psh)o = (Psh)i, and 
the condition of equal velocity around the CD yield a quartic 
equation for the lab-frame BLF Fsh: 

Y)x^ + 2(9 - l)g{Y - g^)x' + 2(2 - j)g^{Y - l)x^ 



71 



where: 



+2(7 - l)g{g^Y - l)x + 9(1 - g^Y) = 0, 



(11) 



9 ' 



9 



<x < g 



(12) 



where Fo and F^ are the initial BLFs of the two shells and Y 
is the ratio of the pre-shocks densities of the slow and fast 
shells. Eq.(jl^) is solved numerically. In the range of interest 
and for any physical choice of parameters Eq.(|ll|) yields only 
one solution. 



• The shell widths after the shocks. The compression ra- 
tio p/peh is obtained by the jump conditions (^ and (jlo|) 
(assuming a value for the adiabatic index 7 = 4/3) and: 



sh,o 



Psh,o 



(13) 



where Ai and Ao are the pre-shock values. 

• The internal energies iJRs and Eps for the RS and FS. 
The equality of the two shocks velocities and of the comoving 
internal energy densities yields (in the lab frame) to: 



-Ers 
Efs 



R. 



(14) 



Thus the energy in the RS is Ens ~ ^/(l + i?) x Ein, and 
in the FS is Eps = 1/(1 + -R) x _Ein, where Ein is the total 
internal energy generated during the collision (eq. 

• Knowing the shell width and the shock velocity we can 
then estimate the shell crossing time of the shock Stcr = 
A/\vsh — vo\ where Vsh and vo are the post-shock and pre- 
shock material velocity, respectively. 



2.4 Shocked fluid parameters 

The parameters determining the synchrotron and inverse 
Compton emission are the magnetic field strength B and 
the electron energy distribution in the shocked fluid. 

Given the uncertainty in determining these quantities 
directly from the physics of relativistic shocks, these are 
parametrized by dimensionless parameters deflned as: ea 
and £0, which measure the fraction of the co-moving in- 
ternal energy density U!^^^ stored in magnetic fleld and 
in electrons Uc, respectively (the apex indicates quantities 
measured in the co-moving frame). 

The co-moving magnetic fleld is then deflned by 
B'V(8^) = eefy,',,: 



B' = — 

Fsh 



2£BgBh 
i?2A ' 



(15) 



and is assumed to be randomly oriented in space. 

Following the results on non-relativistic shocks, and the 
recent results on relativistic ones (Bednarz & Ostrowski 
1998; Kirk et al. 2000), we also postulate that the elec- 
trons are accelerated (instantaneously) behind the shock as a 
power law A''(7) oc j^^ with 7 > 7b, and p > 2. For an equal 
number of protons and electrons, the number density of non- 
thermal electrons is: nj, = CcPsh/mp, where is the fraction 
of electrons effectively accelerated. Since — ecUsh, the 
minimum RLF 7b corresponds to: 



7b 



p-2 
p-1 



1 + ^EkilUx 



1837^ ^F,b. 

P - 1 Cc 



(16) 



In conclusion, by studying the wind dynamics and shock 
hydrodynamics we can determine for each collision the en- 
ergy density in non-thermal electrons and in magnetic fleld 
and estimate the minimum RLFs 7b of the electron energy 
distribution injected by the shocks. 



Once Fsh is known we can evaluate: 
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3 THE SPECTRUM 

Having estimated such quantities, we can therefore calculate 
the spectrum of the radiation emitted during each collision. 
In order to achieve that wc adopt some simplifications: 

• The emitting zone is assumed to be homogeneous, and 
the plasma is embedded in a tangled magnetic field. Its co- 
moving volume is V = TT^{P R'^ AR' , where tp is the half- 
opening angle of the jet, assumed conical. 

• The relativistic particles are treated as having the same 
energy distribution throughout the entire emitting region. 
In a more realistic approach, the particles would be acceler- 
ated at the shock front, where most of the energetic electrons 
would bo located, while cooled particles would be progres- 
sively more distant from the shock front. However, as we are 
interested in the "average" spectrum, we assume to collect 
radiation emitted by the entire shell. Clearly such simplifi- 
cation docs not allow to consider changes in spectral details 
occurring on a timescale faster than the light crossing time 
of a single shell. 

• There is a source of soft photons external to the jet. We 
identify it with the emission reprocessed in the broad line 
region (BLR). Therefore wo consider a luminosity Lcxt = 
aZ/disk is produced within /Jblr, corresponding to a radi- 
ation energy density (as measured in the frame comoving 
with the shell) given by (e.g. Ghisellini & Madau 1996): 



17 aLdiskr^ 



12 47r7?|LRC 



(17) 



For simplicity, this seed photon component is considered to 
abruptly vanish beyond i?BLR. 

• We assume that the particle distribution results from 
the continuous injection and cooling processes, where par- 
ticles are injected with a power law energy distribution of 
index s = {p — 1) between 7b and 7max (as measured in 
the comoving frame). The corresponding power injected in 
the form of relativistic particles is Le. The equilibrium par- 
ticle distribution A'^(7) is then found through the following 
procedure. 

We calculate the (comoving) random RLF 7cooi of those par- 
ticles with a cooling time comparable with the shock cross- 
ing time Across and above 7cooi the distribution is taken as a 
power law of index p. Below 7cooi there are two cases: 

i) if 7cooi > 7b, we assume that A''(7) oc 7~^ between 7b 
and 7cooi, and A''(7) = for 7 < 7b; 

ii) if 7cooi < 7b, -^^(7) oc 7~^ between 7cooi and 7b, and 
N{'y) = for 7 < 7cooi- 

The normalization of N{'y) is determined according to 
whether particles with 7 ~ 7b can or cannot cool in the 
timescale taoss'- 

a) if most of the power injected in relativistic particles is 
radiated in a short timescale (< f cross, this is equivalent to 
demand 7b > 7cooi) we apply a luminosity balance condition: 



Vmec' / jN{'y)d'y, 



(18) 



where V is the emitting volume and 7 is the total radia- 
tive cooling rate, which includes synchrotron, synchrotron 
self Compton (SSC) and external Compton (EC) losses. In 
order to allow for the possibility of having multiple Comp- 
ton scatterings, especially when 7b is not very large, we first 
calculate how many scattering orders nic can take place be- 



fore the Klein-Nishina limit is reached (we consider up to 
10 orders of scatterings). Taking 7b as the relevant energy, 
the number of scatering orders in the Thomson regime are 

ln(7b/a;B) 



nic = 



- 1 



(19) 



ln(47^/3) 

where xb = hi/B/{me(?) and 1/3 = eB/i^-Kmec) is the (non- 
relativistic) Laxmor frequency. We define Ue = Le/{i^R^c) 
and introduce the Comptonization parameter y as: 



\(jtAR 
4 



7 



7 A^(7)d7. 



(20) 



Here A7i' is the thickness of the shell at the time of the 
collision. For semplicitly, we neglect its further expansion 
while it is radiating. With this definition of y, Eq. (18) can 
be re-written as: 



t^ext I 2 . 

l + — +y + y + 



...y 



b) if most of the power is not radiated in t, 
retain their energy and 



7 



rrieC V / N{-y)-yd'y. 



(21) 

cross, the electrons 
(22) 



Assuming the same broken power-law shape for A'^(7) we 
can then determine its normalization. 



3.1 Synchrotron emission 

To calculate the synchrotron emission we adopt the pitch 
angle averaged synchrotron emissivity for a single electron 
given in, e.g., Ghisellini & Svensson (1991): 

, 81/3 (TtcUb ^2 

TT Vb 

[K,/,{t)K,„{t) - I [Kl,S) - Kl/S)] } ■ (23) 

Here t = i^/(37^i^b) and K\, is the modified Bessel function 
of order h. The total monochromatic synchrotron luminosity 
(in the comoving frame) is then found by integrating j(7, v) 
over the particle distribution and the source volume. 



3.2 Synchrotron self— absorption 

The synchrotron self absorption spectrum produced by a 

power law energy particle distribution with a low energy 
cut-off is F{v) oc instead of the familiar u''^'^ . The type of 
distribution can result from the injection of particles above 
some minimum energy and slow cooling (with respect to 
Across). This case is therefore relevant to us. Here we treat 
synchrotron self-absorption exactly, making use of the con- 
cept of "synchrotron cross section" 0-3(7,1') given in Eq. 
(2.17) of Ghisellini & Svensson (1991), which we report for 
completeness: 



o-s(7, 



8V37r^ 
15 



-—\K 
B 76 L-^- 



4/3 



(t)-/^?/3(t)]. 



(24) 



where e is the electron charge and a tangled magnetic field 
and an isotropic distribution of pitch angles are assumed. 
The synchrotron absorption optical depth is then given by 



Tt = AR' / a,i^,y)N{y)d-y. 



(25) 
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Note that in this formalism the optical depth can be deter- 
mined also when A''(7) has abrupt cut-offs. 

The resulting monochromatic synchrotron luminosity is 
isotropic in the comoving frame and is calculated according 
to the usual radiation transfer equation: 



j(7,i/)iV(7)d7 



1 



(26) 



3.3 Inverse Compton emission 

We calculate the inverse Compton emission if it occurs in 
the Thomson regime, while neglect the Klein-Nishina regime 
by adopting a scattering cross section equal to the Thom- 
son one for < 3/4 and zero otherwise, x — hv / (rricC?) 
is the dimensionless energy of the incoming photon. The 
5-function approximation for the single particle spectrum 
is adopted, i.e. it is assumed that the particle emission is 
peaked at Xc — (4/3)7^a;. With these simplifications the re- 
sulting emitted monochromatic power is 



Lc{Xc) 



aTcVrricC 



Ni^)^n{x)dx; 7 = (^) (27) 



where n{x) is the density of seed photons. Eq. ( |27[ ) is cal- 
culated for each order of Compton scattering. To calculate 
Lc,n{xc) of the Compton order n we must then use the pho- 
ton density produced by the previous order, n(a::)n-i. For 
the first order, the seed photons are the sum of the locally 
produced synchrotron and the externally produced photons. 
The spectrum of the external seed photons (in the comov- 
ing frame) is assumed to have a blackbody shape peaking at 
frequency xq = Fsiinc, where we always assume Pnnc ~ 
Hz. 

As high energy 7-rays can interact with softer photons 
producing electron-positron pairs, we take into account the 
absorption effects due to this process. However in general 
this is found to be unimportant because no large compacte- 
nesses are involved. 



3.4 Beaming 

Finally let us consider the effect of relativistic beaming on 
the observed radiation. The emitting plasma is moving at a 
BLF F, intermediate between the original BLF of the two 
colliding shells, computed according to the details given in 
Section 2. The radiation produced by each colliding shell is 
then received by the observer as 

^obs^^obs-j = 5^L{u°''yS), (28) 

where 5 = [F— (F^ — 1)^^^ cos6]^^ is the Doppler (or beam- 
ing) factor and 6 is the viewing angle (with respect to the jet 
axis). The luminosity thus calculated corresponds to what 
the observer would estimate assuming isotropic emission. 



4 NUMERICAL SIMULATIONS 

In this section we analyze the general properties of the 
macro- and micro-physics of the plasma during the jet evo- 
lution. Numerical simulations are mandatory in order to test 



the model, given the difficulty in understanding and disen- 
tangle the effect of each parameter on the global properties 
of the flow and the produced radiation. 

The number of parameters that sets the kinematical 
and radiative evolution of the jet flow is relatively small. 
More precisely, the range of BLFs (Fm -r- Fm), the average 
kinetic luminosity of the inner engine Lw and the time vari- 
ability timescale tv deflne the kinematical evolution, while 
the equipartition parameters ec and eb and the fraction of 
accelerated electrons control the radiative properties of 
the plasma. The former ones are reasonably constrained, 
since 10 F ^ 25 and Lw ~ 10** erg s"^ are observation- 
ally inferred and the variability time can be estimated as 
tv ~ -Ro/c ~ 10* s. On the other hand, it is more difficult 
to quantify the latter ones and they can be basically con- 
sidered as free parameters. In the simulations presented in 
this paper, they have been specifically adjusted in order to 
reproduce the general observed behavior of the OVV source 
3C 279, as: = 0.5, eb = 0.004 and Co = 0.04. 

Finally, the external radiation field (see above) is con- 
strained by the radius of the BLR, -Rblr = 5 x lO^^cm, and 
by the luminosity of the disk and corona radiative compo- 
nents assumed to be L = 10*^ erg s"^. 

In order to illustrate the difficulty in predicting the re- 
sults of changes in the above parameters, let us consider, as 
an example, the effect of variations (decrease) in the lower 
BLF Fm. On one side as the average speed of the shells is 
smaller, they must be more massive to maintain the same 
wind luminosity. This implies that the energy per baryon 
generated in each collision is lower and thus a smaller 7b 
is expected. On the contrary, the dynamical efficiency in- 
creases thanks to the larger ratios of the BLFs of the in- 
volved shells, leading in turn to larger internal energy, in di- 
rect competition with the former effect. Moreover, a smaller 
Fm causes the collisions to start earlier, when the shell vol- 
ume is smaller, with the consequent creation/amplification 
of a larger magnetic field. In summary, even a limited change 
in a single parameter causes several (often competing) effects 
in the kinematical and radiative evolution of the flow. There- 
fore in the following we will examine such consequences 
through numerical simulations of the flow evolution. 



4.1 Output parameters 

A first key quantity calculated through the simulations is the 
distribution of collision radii: this in fact controls the aver- 
age value of the magnetic field strength (which depends on 
the shell volume, i.e. on the shock distance), and the num- 
ber of collisions occurring inside the BLR, which in turns 
determines the relative importance of the radio-optical to 
7-ray flux. 

Figure |^ shows the distribution of radii at which colli- 
sions take place. The limit at small distances directly follows 
from the choice of the injection parameters, and its typical 
value can be approximated as: 



Rn 



2al 



(a- 



7 X 10'" cm. 



ar 



Fm 

fT 



(29) 



for an average time between two ejections tv = 10 s, Fm = 
25 and Fm ~ 10. The merged shells collide again and again 
at larger distances up to ^ 10^" cm, where the differences 
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log(R) cm 



Figure 2. The histograms show the distributions of radii of the 
collisions between shells that have; never collided before (i.e. with 
no internal energy), only one has collided before, both have col- 
lided before. 



among the shell velocities are completely smoothed out and 
the wind can be considered uniform. 

The overall duration of the simulation (tw ~ 10*s, corre- 
sponding to 3x10* shells) has been chosen in order to reach 
a virtually steady state for the wind, i.e. with collisions tak- 
ing place at all radii smaller than ~ 10^" cm. As shown 
in Fig. ^ in the inner part {R < 10^* cm) the collisions 
are mostly between shells which have never collided before, 
while shells colliding at larger radii (10^*-10^^cm) are the 
result of previous mergings. About 50 per cent of the shocks 
occur inside the BLR (for Rblr = 5 x 10^^ cm), while only 
7.2 per cent take place at radii R > 10^® cm. 

The net amount of internal energy Ein released in each 
interaction depends on the dynamical efflciency and on the 
total kinetic energy of the shells available before the colli- 
sion. As illustrated in the upper panel of Figure ^ the shell 
BLFs tend to average out during the wind expansion, lead- 
ing to a strong attenuation of the dynamical efficiency: this 
can be as high as 10 per cent for the inner collisions but 
decreases to less then 0.1 per cent at the largest radii. How- 
ever, the outer shocks involve larger kinetic energies, as the 
shells are much more massive - being the result of many pre- 
vious collisions. As a consequence, the average net amount 
of internal energy generated by each interaction is roughly 
constant with radius, although as few shells are left in the 
large scale wind, their rate of collisions is low. 

As discussed in the previous section, the internal energy 
Ein is shared among protons, electrons and magnetic field. 
The central panel of Fig. ^ shows the value of the minimum 
RLF of the accelerated electrons, 7b, while the lower panel of 
Fig. ^ represents the magnetic field strength B for the same 
collisions. In some cases, the efficiency of the collision is so 
small that the minimum Lorentz factor of the accelerated 
electrons is lower than 2. Also the magnetic field is usually 
small in these cases. These collisions are not considered in 



iog(Yj 



log(B) 




17 18 19 

log(collision radius) cm 



20 



Figure 3. Evolution of the post-shock shell parameters as a func- 
tion of the collision radius. The upper panel shows the BLF of 
the merged shell after each collision; the central panel represents 
the minimum RLF of the electrons 7^, (for te = 0.5; = 0.04); 
the lower panel shows the value of the generated magnetic field 
(eb = 0.004). FS and RS are considered independently. Only one 
tenth of the points is displayed for clarity. Circles refer to colli- 
sions in which none of the shells has ever collided before, pluses 
to collisions in which one of the shells has collided, while triangles 
indicate collisions between shells that have both collided before. 



the radiative part of the code, since the cooling time of the 
electrons is so long that adiabatic losses are overwhelming. 
In any case, the minimum Lorentz factor of the electrons 
always excedes the thermal motion of the fraction of non 
accelerated ones. Since Ein is approximately constant with 
radius, the value of B strongly anti-correlates with the col- 
lision distance approximately as B cx R~^'^, being inversely 
proportional to the square root of the volume. 

We can calculate what fraction of the energy produced 
at each location of the jet is used to accelerate electrons 
to relativistic energy. In the absence of non-radiative losses, 
this energy will be eventually entirely radiated. We then call 
this quantity, which is independent of the specific radiation 
mechanism and the timescales involved, the "differential ra- 
diative efficiency". The logarithmic differential efficiency, i.e. 
the efhciency per unit logarithmic radius, is shown in Fig. ^ 
where we can see that the most radiatively efficient region 
is the inner jet (between 10^®'^ and 10^* cm), where the 
collision rate is the highest. 

The corresponding integrated efficiency, i.e. the fraction 
of Lw radiated from the beginning of the jet up to a given 
radius, is shown in Fig. ^: it saturates at i? ~ 10^® cm at a 
value of ~ 6 per cent. Such value can in principle depend on 
the input parameters, but an efficiency of a few per cent is 
a robust result unl ess a very broad dis tribution of the initial 
BLFs is involved (Beloborodov 2000). Note that the finite 
duration of the computation does not affect the efficiency 
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Figure 4. The radiative efficiency versus the collision radius. The 
solid line refers to the global efficiency, i.e. the fraction of the to- 
tal kinetic energy of the wind {E^ = Lw x t^) radiated on scales 
smaller than a given radius; the shaded hystogram shows instead 
the differential efficiency, namely the fraction of -Ew radiated for 
a given radius interval. The hystogram amplitude has been mul- 
tiplied by a factor of 10 for clarity. The vertical line indicates the 
radius assumed for the extension of the BLR. The cone-like insert 
in the upper part of the figure shows a grey-tone representation 
of the differential efficiency of the jet. The darker the color the 
higher the efficiency (cfr. the shaded hystogram). 



estimate, since the difFerential value approaches zero at radii 
~ 10^^ cm, much smaller than the maximum radius reached 
in the simulation. 

From the micro-physical properties of the plasma 
shown in Fig. |^ it is also possible to compute the typical 
frequencies of the emitted photons. For the inner collisions 
{R < -Rblr) synchrotron, SSC and EC are responsible for 
the observed radiation. At larger radii, external to the pho- 
ton bath provided by the BLR re-processing, the first two 
mechanisms dominate the emission. The evolution of the 
peak frequencies corresponding to the three radiative pro- 
cesses along the jet is shown in Fig. |^, where their values at 
each radius have been obtained by averaging over all of the 
collisions occurring at that radius. The synchrotron peak fre- 
quency I'syn cx 7^5 decreases from ~ 10^'' Hz for the initial 
collisions to ~ 1 GHz at pc scale. The SSC peak frequency, 
given by ussc oc j'^Vsyn, varies from ~ 10^° to ~ 10^^ Hz, for 
R increasing from 3x10" to 10^** cm. Finally the EC typical 
frequency uec ~ ^mlhi^inc ~ 10^^ Hz (for i^unc ~ 10^^ Hz). 

Let us now consider the evolution of the effective radia- 
tive efficiency, which is estimated by comparing the radiative 
time scales St^ of relativistic electrons (corresponding to the 
cooling due to the three above processes) with the time a 
shock takes to cross the shell width (i.e. the energy injec- 
tion time) Stci, and the angular spread time, Stg (see Fig. |^. 
StcT increases from a few hr at 5 x 10^® cm to ~ a week at 



10^* cm, because of both the increase in the shell width and 
the slowing down of the shock velocity during the expansion. 
The angular spread time refiects the different distances trav- 



eled by photons propagating along different directions. Since 
the observed radiation is beamed in a solid angle ~ l/F, Ste 
corresponds to the delay in the arrival time of photons emit- 
ted at an angle 1 /F with respect to those propagating along 
the line of sight: Stg « R/2cr^. This varies from a few hr to 
~ a week, analogously to the shell crossing time-scale. The 
cooling time 5t^ is determined by the most efficient process 
among synchrotron, SSC and EC. Inside the BLR the latter 
dominates, with a cooling time of less than a hr, over the ~ 1 
day time scale of synchrotron and SSC emission. Outside the 
BLR, the synchrotron process dominates the emission, with 
strongly radial dependent cooling times tsyn ~ 1 yr Note 
that inside the BLR the cooling is so efficient that electrons 
emit most of their energy immediately after being acceler- 
ated and the pulse duration is due to a combination of the 
shock crossing time and geometrical effects. At these radii, 
the radiation is mostly emitted (via EC) at high frequencies. 



10^ 



10 Hz. Outside the BLR instead electrons cool on 



a time scale much longer than Sta, and thus the pulse du- 
ration is set by the synchrotron cooling time and the peak 
frequency is in the mm range. 



5 RESULTS 

Via the numerical simulations the full time-resolved spectral 
behavior of the source can be determined. In this section, we 
directly perform comparisons of the model output quantities 
with observational data relative to 3C 279, by decompos- 
ing the information into (snapshots or time averaged) spec- 
tra and light curves. The full time-dependent behavior of 
the simula ted source can be also examined through anima - 
tions (see 



http:/ / www. merate. mi. astro, it/^ lazzati/ 3 C279/ ) 



in which the temporal and spectral evolutions arc simulta- 
neously shown. 



* Clearly the synchrotron and SSC cooling times increase with 
radius, since 7^, and B decrease during the expansion. 
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Figure 6. Time scales determining the duration of the observed 
pulses versus the collision radius. The dashed line refers to the 
angular spread time, the solid line to the time a shock takes to 
cross the shells, and the dotted line to the electron cooling time. 
The circles represent the pulse duration (for clarity we plot only 
one tenth of the efficient collisions). 
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The first step of the simulations (see the previous sec- 
tion) determines a set of physical parameters for every col- 
lision. In a second step these are then used to produce the 
Spectral Energy Distributions (SED) emitted in each colli- 
sion (one SED for the forward and an independent SED for 
the reverse shock) as described in Section ^. Different spectra 
are generated in different collisions, but general behaviors 
which describe the overall shape of the SED as a function of 
radius can be found. As already stressed the most important 
condition that affects the SED shape is whether the collision 
occurs inside or outside the BLR. In Figure |^ we show the 
averaged spectra of collisions taking place within a given 
range of radii. Shocks inside the BLR produce spectra dom- 
inated by the EC process, through the scattering of electrons 
with energy 7b meC'^ on BLR photons. The synchrotron peak, 
due to the same electrons, is located in the near infrared. A 
minimum occurs in the optical-UV range, while the X-ray 
band is characterized by a power-law spectrum. Outside the 
BLR (lacking the EC component) the SED becomes more 
and more dominated by the synchrotron process: immedi- 
ately outside the BLR the Compton parameter y ~ unity 
and the spectrum is the sum of the synchrotron and the 
first two Compton orders; as the radius increases, both the 
relativistic particle density and their average (squared) en- 
ergy decrease, leading to progressively smaller values of y (if 
only radiative losses are considered). At the same time the 
synchrotron self-absorption frequency decreases with radius. 



Figure 7. Average spectra of collisions. Solid line show averages 
for collisions happening inside the broad line region _R < 5 X 
10^^ cm, dashed line shows collisions happening at intermediate 
distances 5 X 10^^ < R < 10^® cm and dot— dashed line show 
collisions at large radii {R > 10^^ cm). Upper panel show the 
average spectrum of the shells without taking into account the 
duration of the emission. The Figure show the average spectrum 
of a single shell in the given range of radii. In the lower panel, 
instead, the duration of the emission of each collision has been 
taken into account. The plot shows the average instant spectrum 
of the whole section of the jet. 

sum of all the spectra produced by shocks simultaneously 
active in the observer frame, taking into account the radius 
at which each collision is occurring and the travel path of 
photons to reach the observer. 

We therefore associate to each spectrum a photon pulse, 
with a start time (the time at which the collision begins) 
and a duration (set by the combination of the timescales 
discussed in Section Uf). The shape of such pulse is trian- 
gular if its duration is determined by geometrical effects 
and/or by the shock crossing time (with a FWHM equal 
to the timescale shown in Fig. these are generally asso- 
ciated to collisions inside the BLR. If the pulse duration is 
instead set by cooling the pulse rises linearly, as in the pre- 
vious case, but after the maximum it decays exponentially 
[oc exp(t/3tcooi(7b))]: in this case the rise and decay times 
can be very different. 



5.1 The observed spectrum 

We have so far examined the role of the various emission 
mechanisms at different distances along the jet. The ob- 
served spectrum is the convolution of the emission of all pho- 
tons arriving at the observer simultaneously: since these are 
produced at different distances, it is necessary to take into 
account the different light propagations times when sum- 
ming them up. In other words, the resulting spectrum is the 



5.2 The spectrum of 3C 279 

Examples of the resulting spectra are compared to data from 
three different simultaneous multiwavelength campaigns of 
3C279 in Figure ^. Note the remarkable variability, espe- 
cially at 7-ray energies. 

Within the model examined here the source high state 
corresponds to a time interval in which at least one of the 
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Figure 8. Example of the spectra predicted by our model com 
pared with the SED of 3C 279 corresponding to different obser 
vational campaigns. 

Simulated lightcurves 



10 " r 



10 



10 



10 



10 



10 



g 10" 



M 10 



= 10 



-10 



-11 



12 



-10 



-11 



12 



13 



-10 





10 



10 



10 



10 



-11 



12 



-10 



-11 



Gamma Rays (400 MeV) 



X-Rays (1 KeV) 



Optical (V) 




Infrared (1 mm) 



100 



300 
Time (days) 



300 



Figure 9. Light curves at different frequencies. The 7-rays, X- 
rays and Optical Ught curves vary on different time scale, with a 
minimum value of few hr, and the IR one varies on few months 
time scale. 
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Figure 10. Results of cross correlations between 7-rays and X- 
ray (solid grey line), optical (solid dark line) and infrared (dashed 
line) light curves. The X-rays and optical emission vary simulta- 
neously with the 7-rays, while the IR follows the 7-rays emission 
with a delay of a ~ 40 days (the vertical dotted line). 



active regions is inside the BLR. The lower line shows the 
SED in a quiescent state of the source, i.e. when all of the 
active collisions are located outside the BLR. The differ- 
ence between the two spectra manifests itself mainly in the 
high frequency range, where the internal collisions, that have 
a shorter pulse duration, mostly contribute. Although the 
number of external collisions is almost equal to the number 
of the inner ones, the former produce longer pulses: it follows 
that the number of those simultaneously active is higher and 
hence the statistical variability is suppressed. 

A different way to examine the results is to determine 
light curves in given energy bands: Figure ^ shows the flux 
temporal behavior in the 7-ray, X-ray, optical and millime- 
ter ranges. The high frequency fluxes (down to optical en- 
ergies) are highly variable (again since they are produced 
mainly by the inner collisions), while the millimeter flux 
comes from outer shocks and thus the corresponding hght 
curve is much smoother. 

In order to study the presence of delays in the variabil- 
ity at different wavelengths, we added "secular" variations 
in the kinetic luminosity injected by the inner engine. In 
particular, this is modulated as a sinusoidal function with 
period ~ 1 year, spanning a factor of five in luminosity (be- 
tween the minimum and the maximum). Figure |l^ illus- 
trates the results of the cross correlations between the 7-ray 
and the other frequencies light curves. Once again the plot 
shows that 7-ray, X-ray and optical photons are produced 
in the same region (being closely correlated, with no time 
delay), while the IR emission lags by ~ 40 days, as the shells 
that produce high frequency photons inside the BLR have to 
travel to larger radii {R ~ 10^^ cm) to produce IR emission. 
The expected observed delay is given by: 

AR 



= — ~ 38.5 Ai?i9 rr days. 



(30) 
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which indeed well describes the results of the simulations. 

Finally note that time delays of the order of ~ hr might 
occur between 7-ray and X-ray/optical peaks even if they 
are produced in the same region. However, the study of such 
delays would require to take into account the different light 
travel paths within a single region (Chiaberge & Ghisellini 
1999), and a time dependent treatment of the shell emission, 
taking into account the variations in the magnetic field and 
electron distribution during the cooling. Such a treatment is 
beyond the scope of the present work and will be examined 
in future developments of the model. 



6 DISCUSSION 

Highly relativistic flows dominate the dynamics and emis- 
sion in both Gamma-Ray Bursts and blazars. It has been 
proposed that internal shocks, produced by the interaction 
of flow components injected with different bulk speeds, can 
be responsible for the observed emission and variability in 
both classes of sources (Rees 1978; Rees & Meszaros 1994). 

Here we have considered in some details this 'internal 
shock' scenario in the context of AGN jets, self-consistently 
accounting for the dynamics and the radiative properties 
of the flow. After describing the kinematic and radiative 
assumptions of the model, we have examined the results 
of numerical simulations which allow us to study the pre- 
dicted energy deposition on different jet scales, the rele- 
vant timescales, the corresponding dominating emission pro- 
cesses, the spectral evolution and predicted variability, the 
correlations among the emission at different frequencies. 

The key interesting features of this scenario are the rel- 
atively low radiative efficiency (of order of a few per cent) 
which well accounts for the dissipation of kinetic energy in 
blazars, as higher dissipation rates would create major en- 
ergetic problems (in order to carry sufficient power up to 
the large scale lobe structures). Furthermore, although the 
radiative dissipation occurs on all jet scales, as indeed ob- 
served, the bulk of it is localized at tenths of a pc, on the 
BLR scale, in agreement with the requirements of fast vari- 
ability and transparency to 7-rays. 

The results of the simulations, based on input parame- 
ters suitable to qualitatively reproduce the observations of 
a powerful radio-loud quasar, such as 3C 279, are indeed 
satisfactory. As expected the plasma in the inner part of the 
jet cools predominantly via EC emission, while on larger 
scales synchrotron cooling dominates. Qualitative agree- 
ment is also obtained in terms of light curves, and corre- 
lations among difi'erent spectral bands are predicted. Also, 
an asymptotic value of the bulk Lorentz factor F ~ 15 is 
achieved at distances ~ pc. 

While this paper has been aimed at presenting the 
model and study its general applicability to blazar jets, 
the final goal is to determine which initial parame- 
ters/conditions are required in order to produce the whole 
of the blazar phenomenology and energetics. In particular, 
the model significantly constrains the typical Lorentz fac- 
tors (and their distribution) achieved through acceleration 
on the inner scales, and the typical timescale of plasma in- 
jection. In this respect, also the determination of the duty 
cycle of the 7-ray emission in blazars is of chief importance, 
as it can be used to constrain the rate and Lorentz factor 



distribution of the injection of energy/plasma in the rela- 
tivistic flow. Note that although more accurate measures of 
7-ray spectra and the detection of a much larger number of 
7-ray loud sources will require the capabilities of GLAST, 
tight constraints on the duty cycles will be available much 
sooner thanks to the launch of AGILE. The model presented 
here will provide (for the flrst time) the possibility of directly 
modelling the time dependent broad band spectral evolution 
of these sources following the (time dependent) energy in- 
jection. The relationship between flares at high energies and 
the birth of superluminal radio knots and/or radio flares will 
be also explored. 

Therefore issues associated with the jet energy content 
and the relative conversion of bulk energy into internal (par- 
ticle and fleld) energy at relativistic shocks can be also con- 
strained by the comparison of the model predictions with 
observations. 

A further aspect we intend to examine is the propaga- 
tion and dissipation of jets on larger scales, where however 
uncertainties are large. On the 3> pc-scales, signiflcant en- 
trainment might occur, possibly leading to the formation of 
external shocks and the development of slower outer lay- 
ers around jets. Another issues is the relative importance of 
creation vs amplification of an existing seed magnetic field: 
if the latter process prevails, it will probably lead to larger 
magnetic fields on the kpc scale jet, and consequently to 
faster synchrotron cooling. Inverse Compton cooling at kpc 
scales is instead enhanced, with respect to what calculated 
here, by the interaction of with the cosmic background radi- 
ation, which is seen relativistically boosted in the rest frame 
of the fast moving parts of the jets (e.g. Celotti, Ghisellini 
& Chiaberge, 2000). 

We finally stress that it is of primary importance, for the 
understanding of the flow injection and its connection with 
the central accreting engine, to determine the key param- 
eters which regulate the observed spectral and luminosity 
trend, from high power, low energy peaked, flat spectrum 
radio quasars, to weak high energy peaked BL Lacs (Fossati 
et al. 1998) . This will be the content of a forthcoming paper. 
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